Fermionized photons in an array of driven dissipative nonlinear cavities 
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We theoretically investigate the optical response of a one-dimensional array of strongly nonlinear 
optical microcavities. When the optical nonlinearity is much larger than both losses and inter-cavity 
tunnel coupling, the non-equilibrium steady state of the system is reminiscent of a strongly correlated 
Tonks-Girardeau gas of impenetrable bosons. Signatures of strong correlations are identified in the 
absorption spectrum of the system, as well as in the intensity correlations of the emitted light. 
Possible experimental implementations in state-of-the-art solid-state devices are discussed. 
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Strong correlations in quantum many body systems 
give rise to a number of striking phenomena and states 
of matter, ranging from superfluidity of liquid Helium 
and superconductivity of metals to fractional quantum 
Hall effect in two-dimensional electron gases. Most of 
the work in condensed matter physics so far focused on 
systems close to thermodynamic equilibrium. It is known 
however that even richer range of behaviors and phases is 
exhibited by systems where the steady state arises from a 
dynamical balance between drive and dissipation rather 
than from a thermodynamic equilibrium condition [l|, Q • 
The study of non-equilibrium phenomena has been his- 
torically confined to classical models, such as driven dif- 
fusive lattice gases |3|] and models of hydrodynamic sta- 
bility ■ Another well-studied, prototypical example of 
a non-equilibrium phenomenon is lasing Q. In contrast, 
non-equilibrium phases of quantum many-body systems 
remained largely unexplored. Early work in this direc- 
tion has addressed dissipative magnetic chains [6|[, driven 
electron gases , and nonlinear optical systems [1, [^, ■ 

Several papers have recently investigated the possibil- 
ity of observing strongly correlated quantum phases in 
gases of interacting photons [lH , including Mott insula- 
tor to superfluid transition in cavity arrays I4I, inter- 
acting spin models Il3l|. the Tonks-Girardeau gas in a 
waveguide geometry (l4|. and quantum Hall states [l5| . 
The focus of almost all this prior work has been the real- 
ization of well-known and thoroughly studied equilibrium 
many-body systems in photonic media. 

In this Letter we report the theoretical investigation 
of a genuine non-equilibrium state of strongly correlated 
photons. We study the spectroscopic signatures of strong 
interactions in an array of cavities mutually coupled by 
tunneling and driven by a coherent laser field. The prop- 
erties of the non-equilibrium steady state are studied as a 
function of the pump frequency and intensity for a range 
of system parameters. Specific attention is devoted to 
the strongly nonlinear case, where non-trivial quantum 



correlations appear between photons indicating the on- 
set of a Tonks-Girardeau gas of "fermionized photons" . 
Most significantly, because the quantum correlations in- 
side the many-cavity system directly transfer to the emit- 
ted radiation [l6| . a much wider range of observables is 
experimentally accessible than in analogous systems of 
ultracold atoms [itI]. 

The system we consider is sketched in the left panel 
of FiglT] Assuming that the dynamics is restricted to a 
single photon mode per cavity, the Hamiltonian can be 
written in the following generalized Bose- Hubbard form: 

H = ^ hujQclci + hU c\c\ciCi — ^ hJ c\c-j-\- 



+ Y^[F,{t)cl+F:{t)c,] , (1) 



where the Ci (cJ) operators destroy (create) a photon 
in the i-th cavity located at x^. indicates next- 

neighbor cavities, C/ > is the (repulsive) single-photon 
Kerr interaction within each cavity, J is the inter-cavity 
hopping energy coming from the overlap of the nearest- 
neighbors cavity fields, and coq is the bare frequency of the 
isolated cavities. The term proportional to Fi{t) accounts 
for the coherent drive by the pump laser: in the follow- 
ing we will restrict to the case of a monochromatic pump 
beam of frequency Wp, wavevector kp, and amplitude Fp, 
i.e. Fi{t) = Fp exp[i(kpXi — ujpt)]. Photons are assumed 
to be lost from the system at a spatially uniform rate 7. 
The master equation describing the time evolution of the 
density matrix p(t) then has the standard Liouville form: 
dtp = -ilH, p] + ^ I]i(2ciPc|' - clc.p - pcjci). 

Let us start by discussing the eigenstates of the N 
boson problem in the absence of pumping and dissipa- 
tion. In the linear regime U — 0, photons occupy single- 
particle states of the hopping Hamiltonian, with an en- 
ergy dispersion e{k) ^ ujq — 2J cos(fc). Wavevector is de- 
fined here as a dimensionlcss quantity: the first Brillouin 
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zone (FBZ) then corresponds to the interval k S [— tt, tt]. 

In the opposite Umit of impenetrable bosons U/ J = oo, 
a generic bosonic iV-body wavefunction . . can 
be exactly mapped onto a fermionic wavefunction by the 
transformation llSl: 



, ZAr)e(cr) . 



(2) 



Here, ii,...,iN are the positions of the N particles, 
a is the permutation that sorts the spatial coordinates 
ii,...,iN into ascending order, and e{a) is the sign of 
the permutation a. This sign guarantees that for any 
wavefunction ip symmetric under the exchange of any 
two particles, the corresponding tpp is anti-symmetric as 
required by Fermi statistics. 
As shown in 
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the eigenstates of the impenetra- 
ble boson problem are in a one-to-one correspondence 
with those of the non-interacting fermionic system, which 
are in turn simply classified by the occupation numbers 
of single-particle orbitals. In the following, we shall 
use the shorthand notation \qi . . .qn) to indicate the 
bosonic eigenstate corresponding to a Fermi wavefunc- 
tion with one particle in each of the gi, . . . , gjv orbitals. 
The (pseudo-) momenta qa=i...N are to be chosen within 
the FBZ, i.e. Qa £ [— 7r,7r]. Both the energy and the to- 
tal momentum of the bosonic state are equal to the ones 
of the corresponding fermionic one, say E = e(<Za) 
and P = qa (as typical of a lattice, momentum is 
here defined only modulo 2tt). On the other hand, the 
momentum distribution of the bosons is not preserved 
by the Bose-Fermi mapping. In particular, the pseudo- 
wavevectors qa of the a — I . . . N fermionic orbitals do 
not have a direct meaning in terms of physical observ- 
ables of the bosonic system and in particular do not cor- 
respond to their physical momentum k (isl . 

In what follows, we shall focus our attention on a neck- 
lace of M cavities with periodic boundary conditions, i.e. 

■ ■ ,ia = 0, . . .) — ■ ■ ,ia — M, . . .), Va. A remark- 
able feature of the Bose-Fermi mapping is that the pe- 
riodicity condition on the iV-body bosonic wavefunction 
does not directly transfer to the single-particle orbitals of 
the fermionized wavefunction: depending on the number 
TV of particles in the system, the fermionic orbitals have 
to fulfill either periodic (if TV is odd) or anti-periodic (if 
N is even) boundary conditions 2^ 21 \ ■ This reflects on 
the quantization of the pseudo- wavevector q = 27rn/M (if 
N is odd) or g = 2'K{n + l/2)/M (if TV is even), where n is 
an integer number. This peculiar quantization rule leads 
to the following explicit form of bosonic wavefunction of 
the lowest N — 2 state \q, —q) with q = n/M: 
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(3) 



In contrast to equilibrium systems, pump and losses in- 
duce transitions between states with a different number 
of photons. This enables one to extract detailed informa- 
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FIG. 1: Left panel: sketch of the system under consideration. 
Right panel: absorption spectra of the mean number n of 
photons present in each cavity as a function of the frequency 
Atjp = tJp — ojo of the pump beam, taken at normal incidence 
fcp = 0. System of M = 5 cavities in the impenetrable boson 
limit U I J = oQ with J/7 = 20. Different curves correspond to 
increasing values of the pump amplitude F^j'^ — 0.1, 0.3, 1, 
2, 3. The vertical dotted lines indicate the spectral positions 
of the peaks predicted by the fermionization procedure with 
different number of photons. The kets |gi . . . (?jv) indicate the 
pseudo-momenta . . . gjv of the occupied fermionic orbitals. 



tion on the microscopic physics of the strongly interact- 
ing photon gas from the spectra of observable quantities. 
These can be numerically calculated by directly finding 
the stationary state of the master equation or by integrat- 
ing it for long enough times via the Monte Carlo Wave 
Function technique [22| . Examples of absorption spectra 
arc shown in Figs[T]and[2Ka,b), where the mean number 
of photons per site n — (clci) is plotted as a function of 
the pump frequency LUp for a pump at normal incidence 
kp — and a weak loss rate 7 <C J,U. 

For very weak driving \Fp\ <C 7, the dynamics is mostly 
restricted to the vacuum state and the = 1, |g = 0) 
state: the resonant driving of this transition is respon- 
sible for the main peak that is visible in all spectra at 
Aujp = — 2J. For stronger driving amplitudes, higher- 



lying A > 1 states start to appear. Photon number con- 
servation implies that states containing A^ photons in the 
many-cavity system is reached by repeated absorption of 
A^ photons from the coherent drive. A generic many- 
body state I/) of energy Ef will then appear in the spec- 
tra as a narrow resonance peak at frequency Wp = Ef/N. 
In the impenetrable boson limit shown in Fig[Tl the po- 
sition of the peaks can be successfully compared to the 
analytical predictions of the Bose-Fermi mapping indi- 
cated by the vertical lines: each peak is associated to a 
set qa=i...N of pseudo-momenta that is compatible with 
momentum conservation P = q^, = 0. Further con- 
firmation of the peak assignments has been numerically 
obtained by looking at their power dependence: the in- 
tensity of a A^ particle peak starts in fact as On 
the other hand, even though the pseudo-wavevectors qa 
completely identify a quantum state, they do not corre- 
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spond to the physical momenta of the bosonic particles, 
so that no direct information can be inferred on the mi- 
croscopic nature of the state from the momentum distri- 
bution of the corresponding emission. 
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FIG. 2: Panels (a,b) Spectra of the population in the k — 
and k = 2tt/3 bosonic modes as a function of pump frequency 
for a fixed pump amplitude -Fp/7 = 0.5 and different values of 
the nonlinear coupling U/J. Panel (c,d): position of the peak 
and relative occupation n{k — 2n /3)/n{k — 0) at the peak 
position as a function of U/J. Blue dashed lines: asymptotic 
values in the weak U/J and strong U/J ^ 1 interaction 
limits. System of AJ — 3 cavities with J/7 — 20. Pump at 
normal incidence kp = 0. 

In order to fully characterize the transition from the 
weakly to the strongly interacting regime, we have per- 
formed systematic numerical calculations for increasing 
values of the nonlinear coupling U/J (Fig^. These cal- 
culations have systematically and in full detail explored 
the M = 3 case and we have checked that our find- 
ings extend in a straightforward way to larger systems 
with more sites. Spectra of the population Uk = (i>\hk) 



-ikj- 



Cj/vM) in the k — and k — 27r/3 modes 



are plotted in the panels (a,b) for various values oiU/J. 
Again, the main peak at Aup = — 2 J corresponds to the 
resonance of the one particle \q = 0) state. The addi- 
tional peak that splits from it as U/j is increased corre- 
sponds to a two-particle final state. The U/J dependence 
of its position and intensity (panels (c),(d)) provides in- 
sight into the microscopic nature of the two-particle state 
in the different regimes. 

In the strong interaction regime {U/J ^ 1), the two- 
particle state is well captured by the lowest N ~ 2 
fermionized state \q, —q) with q — 7r/3. This value of the 
pseudo-momentum imposed by the anti-periodic bound- 
ary condition is a clear signature of the fermionization 
effect and directly reflects into the asymptotic position of 
the peak at Awp/ J = —2 cos(7r/3) = —1. The relative in- 
tensity of the peak on the different modes ?^fc=±27r/3/^^fc=o 
is also in perfect agreement with the analytical predic- 
tions ?ifc=±27r/3/?^fc=o = 1/4. This latter value is obtained 



by discrete Fourier transform of the wave function ([3]). 

In the weakly interacting regime {U/J <C 1), inter- 
actions can be treated within perturbation theory. To 
zeroth order in U/J, the lowest two particle state is a 
factorizable bosonic state with two particles in the k = 
mode and has an energy — 4J. In this limit, it is visible 
only in the fc = spectrum as a peak at AtOp — —2J. At 
the next order, the state energy is blue-shifted by 2U /M 
and the wave function has the following analytical form: 
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The hole around ii — 12 that was complete in the 
fermionized wavefunction of the strong interaction limit 
([3]) is here much less pronounced and its depth scales as 
U / J . Correspondingly, the relative intensity of the peak 
at Aujp ~ —2 J -I- U /M on the k = ±27r/3 momentum 
components grows as rtj.^±27r/3/'T'fc=o — {U/9J)^. Note 
how the interaction-induced blue-shift of the peak even- 
tually saturates for large U/J 3> 1: the kinetic energy 
cost of creating a node at ii = ^2 is well compensated by 
the suppressed interaction energy. 
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FIG. 3: Plot of the auto- (black, solid line) and cross-intensity 
correlations (red, dashed line) of the emission as a function of 
the interaction strength, U/J. Same system as in Fig[2] 

Further information on the microscopic nature of 
the many-body physics of the system can be obtained 
by inspecting the intensity correlations of the emitted 
light [ill quantified by 3(2^(11,12) = {cl^ cj^ Ci., Ci^) / 
Both the auto- {ii = 12) and the cross- {ii ^ 12) cor- 
relations are plotted in FigO as a function of U/ J for 
a pump laser kept exactly on resonance with the {U / J- 
dependent) two-photon transition. In a very weakly in- 
teracting {U/J <^ 1) system, the two-particle peak over- 
laps the one-particle one, and the emitted light inherits 
the Poissonian nature of the pump laser. 

For intermediate values of U/ J ~ 1, the two-particle 
peak is already well separated from the single particle 
peak {U/M ^ 7). The resonant pump laser then selec- 
tively excites the two-particle state and a measurement 
of the intensity correlations of the emission provides de- 
tailed information on the quantum correlations between 
photons in this state. The fact that the system preferen- 
tially contains .A^ = 0, 2 particles rather than 1 is respon- 
sible for the strong bunching. The almost flat shape of 
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the wavefunction (|4|) makes this bunching observable in 
both the auto- and the cross-correlations of the emission. 
On the other hand, when interactions are very strong 
{U / J ^ 1) and the two-particle state has the fermionized 
form ([3]), the cross-correlation remains strongly bunched, 
but the auto-correlation becomes strongly anti-bunched 
as a consequence of the strong on-site interactions. 

The basic building block of a possible device consists 
of a single cavity with a sufficiently strong on-site single- 
photon non-linearity. Although anti-bunched light gen- 
eration has been shown for a number of emitter-c avity 
combinations such as atoms in high-finesse cavities [23j, 
Cooper-pair boxes in super-conducting strip-line cavi- 
ties 27 1, single quantum dots in micro-pillar [2^ and 
photonic crystal defect cavities [2^, scalability to multi- 
cavity devices still remains a challenging task. An impor- 
tant requirement is that the difference between the res- 
onance frequencies of the cavities have to be kept below 
the energy scales U and J. To overcome this difficulty, 
laterally patterned microcavity systems containing quan- 
tum wells as the nonlinear medium appear to be the most 
promising candidates. State-of-the-art micropillar [2^ 
and laterally patterned microcavity [sO] structures are 
expected to provide a sufficiently large polaritonic inter- 
action [i^l for this purpose. In both cases, the emission 
is strongly peaked in the vertical direction and the spec- 
tral position of the polariton mode is to a large degree 
determined by the photonic confinement, which in turn 
can be precisely engineered. Finally, significant tunnel- 
ing between cavities separated by distances UE^to the /Ltm 
range has been theoretically anticipated in [31] for later- 
ally patterned cavities, and also experimentally demon- 
strated in microdisk 84! or photonic crystal [33| coupled 
cavity systems. This makes possible the independent col- 
lection of light emitted from each cavity, in addition to 
the selection in the far field [s^ . 

In conclusion, we have theoretically investigated the 
spectroscopic signatures of a non-equilibrium, strongly 
correlated gas of photons in a one-dimensional array 
of nonlinear cavities with strong photon-photon interac- 
tions and periodic boundary conditions. We have de- 
scribed the imprint of the transition from a weakly in- 
teracting system to a strongly interacting one, on ex- 
perimentally accessible observables. We find that the 
absorption spectra show well-separated many-body reso- 
nances whose position provides an unambiguous signa- 
ture of photonic fermionization; these can be used to 
spectrally isolate a single many-body state. We show how 
the microscopic structure of the non-equilibrium steady 
state of the system can be inferred from the intensity cor- 
relations of the emission. Finally, we believe the present 
work demonstrates the importance of coupled nonlinear 
optical cavity systems in the theoretical and experimen- 
tal investigation of quantum many-body systems out of 
equilibrium. 

IC is grateful to J. Dalibard and Y. Castin for stim- 



ulating discussions at an early stage of the work, 
acknowledge useful discussions with A. Badolato. 



We 



[2] 
[3] 



[4] 

[5] 
[6] 



[7] 



[8] 

[9] 
[10] 

[11] 
[12] 



[13] 

[14] 
[15] 

[16] 
[17] 

[18] 
[19] 
[20] 
[21] 
[22] 

[23] 

[24] 



[25] 
[26] 

[27] 
[28] 

[29] 
[30] 



D. Mukamel, Lecture Notes at NATO school on 
Soft and Fragile systems, 1999. Preprint available as 
cond-mat/0003424 

D. Ruelle, Physics Today 57, 48 (2004). 
B. Schmittman and R. K. P. Zia, Statistical mechanics 
of driven diffusive systems in C. Domb and J. Lebowitz 
eds.. Phase Transitions and Critical Phenomena vol.17, 
Academic Press (1995). 

P.G. Drazin and W.H. Raid, Hydrodynamic Stability, 
Cambridge University Press, Cambridge, 1981. 

H. Haken, Rev. Mod. Phys. 47, 67 (1975). 

P. Werner, K. Volker, M. Troyer, and S. Chakravarty, 
Phys. Rev. Lett. 94, 047201 (2005); T. Prosen and I. 
Pizorn, Phys. Rev. Lett. 101, 105701 (2008). 
M. A. Zudov, R. R. Du, L. N. Pfeiffer, and K. W. West, 
Phys. Rev. Lett. 90, 046807 (2003); R. G. Mani, et al. 
Nature (London) 420, 646 (2002). 

A. Amo et. al., Nature, in press (2008); A. Amo et al., 
preprint larXiv:0812.2748l 

I. Carusotto, C.Ciuti, Phys.Rev.Lett. 93, 166401 (2004). 
D. Gerace, et al, arXiv:0811.3762 (2008). 

M. H. Hartmann, F. G. S. Brandao, and M. B. Plenio, 
Laser & Photon. Rev. 1, 1 (2008). 

M. J. Hartmann, F. G. S. Brandao, M. B. Plenio, Nat. 
Physics 2, 849 (2006); A. D. Greentree, C. Tahan, J. H. 
Cole, L. C. L. HoUenberg, Nat. Physics 2, 856 (2006); D. 
G. Angelakis, M. F. Santos, S. Bose, Phys. Rev. A 76, 
R031805 (2007). 

M. J. Hartmann, F. G. S. L. Brandao and M. B. Plenio, 

Phys. Rev. Lett. 99, 160501 (2007). 

D. E. Chang, et al. Nature Physics 4, 884 (2008). 

J. Cho, D.G. Angelakis, S. Bose, Phys. Rev. Lett. 101, 

246809 (2008).. 

D.Walls, G.Milburn, Quantum optics (Springer, 1994). 

B. Paredes et al., Nature 429, 277 (2004); T. Kinoshita, 
Wenger, and D. S. Weiss, Science 305, 1125 (2004). 
D. Girardeau, J. Math. Phys. 1, 516 (1960). 
Lenard, J. Math. Phys. 5, 930 (1964). 
H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963). 
Carusotto and Y. Castin, New J. Phys. 5, 91 (2003). 



J. Dalibard, Y. Castin, and K. M0lmer, Phys. Rev. Lett. 
68, 580 (1992). 

H. J. Kimble, M. Dagenais, and L. Mandel, Phys. Rev. 
Lett. 39, 691 (1977). 

P. Michler et al. Science 290, 2282 (2000); R. Brouri, 
A. Beveratos, J. P. Poizat, and P. Grangier, Opt Lett. 
25 1294 (2000); K. M. Birnbaum, et al. Nature 436, 87 

(2005) . 

J. P. Reithmaier et al, Nature 432, 197 (2004). 

K. Hennessy, et al. Nature 445, 896 (2007); A. Faraon, 

et al, Nat. Physics 4, 859 (2008). 

A. Wallraff and et al.. Nature 431, 162 (2004). 

A.Verger, C.Ciuti, LCarusotto, Phys.Rev.B 73, 193306 

(2006) . 

D. Bajoni, et al. Phys. Rev. Lett. 100, 047401 (2008). 
O. EI Daif, et al. Appl. Phys. Lett. 88, 061105 (2006); 
D. Lu, et al. Appl. Phys. Lett. 87, 163105 (2005). 



5 



[31] D. Sarchi, 1. Carusotto, M. Wouters, and V. Savona, [33] K. A. Atlasov et al, Opt. Express 16, 16255 (2008). 

Phys. Rev. B 77, 125324 (2008). [34] B. Jassemnejad, A. Bohanrian, J. Lekki, and K. Weiland, 
[32] A. Nakagawa, S. Ishii, and T. Baba, Appl. Phys. Lett. Optical Engineering 47, 053001 (2008). 

86, 041112 (2005). 



